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Abstract 


We report an alternative method to solve second order differential equations which have at most 
four singular points. This method is developed by changing the degrees of the polynomials in 
the basic equation of Nikiforov-Uvarov (NU) method. This is called extended NU method for 
this paper. The eigenvalue solutions of Heun equation and confluent Heun equation are obtained 
via extended NU method. Some quantum mechanical problems such as Coulomb problem on a 
3-sphere, two Coulombically repelling electrons on a sphere and hyperbolic double-well potential 
are investigated by this method. 
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I. INTRODUCTION 


NU method is based on solving hypergeometric second order differential equations with 


special orthogonal functions [l|. Any second order differential equations which have three 
singular points can be converted to a hypergeometric type differential equation by a suitable 
change of variables. In spherical coordinates, for a given potential Schrodinger equation is 
reduced to a generalized hypergeometric-type equation with appropriate coordinate trans¬ 
formation. Then this reduced equation can be solved systematically in order to obtain 
eigenvalues and eigenfunctions. The main equation of NU method is given in the following 


form 


3 ; 


= 0 

a{z) cr (.'^J 


( 1 ) 


where f{z) is a polynomial of at most first-degree, a{z) and a{z) are polynomials of at most 
second-degree and 'ip{z) is a function of hypergeometric-type. The criteria which are related 
to degrees of polynomial coefficients constitute boundary conditions of the method. 

The solution of Schrodinger equation is a fundamental problem in quantum mechanics 
for understanding the physical systems. Since, when the wave function and energy 

eigenvalue are once obtained one has all possible information about the physical properties of 


the system [^. However the solution of the Schrodinger equation is widely complicated for a 
great number of potentials. To obtain an exact solution is only possible for a limited number 
of potentials such as harmonic oscillator potential. Coulomb potential, Kratzer potential, 
etc. j^. The standart solution method is to expand the solution in a power series of 
independent variable and find the recursion relation for all the expansion coefficients ^ . On 
the other hand, the power series solution method involves some detailed calculation techiques 
in order to reach the exact solution. Moreover in many cases usage of approximation methods 
becomes a necessity. However an approximated solution will roughly describe the system. 

In present work the NU method is extended for second order differential equations that 
have at most four singularities. This is carried out by changing the degrees of polynomials 
in the main equation of NU method. In this extended method f{z) will be a polynomial of 
at most second-degree, a{z) and a{z) will be at most third and fourth-degrees respectively. 
Heun equation and confluent Heun equation can be given as examples for the equations 
that provide these boundary conditions of extended NU method. Heun differential equation 


and its confluent forms appear in several branches of physics. For a review see Also a 
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particular solution of Heun equation has been studied in Ref.|6|. 

This paper is organized as follows: 

In second section the NU method is explained in details. 

In third section extended NU method is described and eigenvalue and eigenfunction so¬ 
lutions are given analytically. 

In fourth and hfth sections eigenvalue and eigenfunction solutions of Heun equation and 
confluent Heun equation are derived via the extended NU method. 

In sixth section we report some physical applications about the extended NU method. 


II. NIKIFOROV-UVAROV METHOD 


The foundation of NU method is to solve the second order hypergeometric type differ¬ 
ential equations by means of special orthogonal functions. For a given potential, the radial 
Schrodinger equation is reduced to a generalized hypergeometric differential equation that is 
given by Eq.([T]). Several classes of special functions such as classical orthogonal polynomials, 
spherical harmonics, Bessel and hypergeometric functions are obtained from the solutions 
of this equation [^. By choosing an appropriate (j){z) and taking Eq-dl]) 

can be reduced to a more comprehensible form; 

a{z)y"{z) + T{z)y'{z) + Xy{z) = 0 (2) 


where 

t{z) = f{z) + 2n{z) (3) 


This newly dehned polynomial depends on polynomial 7r{z) which is introduced as; 


a’{z)-T{z) a'{z)-T{z) x , , / x 

7r{z) = -± \ { - — a{z) + ka[z). 


( 4 ) 


In order to determine the polynomial iriz), the parameter k under the square root sign must 
be known and the expresion under the square root sign must be the square of a polynomial 
such that the requirement on 7r{z) namely a polynomial of at most hrst degree has to be 
satished j^. 

By using the property that, all the derivatives of hypergeometric functions are also of the 
hypergeometric type, the solutions of Eq.(l2]) can be generalized. After this generalization 
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the eigenvalue equation can be obtained as follows; 


/ / X n(n — 1 ) ,, , ,, 

X = Xn = -nr (z) - -—-a (z)(n = 0 , 1 , 2 , ...)• ( 5 ) 

For the eigenfunction problem, there is a particular polynomial solution of the form y(z) = 
Uniz) which is given by Rodrigues formula [^; 

r> jn 

where Bn is a normalization constant and p{z) is the weight function and this function 
satisfies the equation {a{z)p{z)y = t{z)p{z). The other function namely (j){z) in the '0(^) = 
4>{z)y{z) transformation is dehned as a logarithmic derivative; 

(j){z) a{z)' 

After determination of all the polynomial parameters, the eigenvalue and the eigenfunction 
solutions of Schrodinger equation can be obtained completely in a systematic way by using 
NU method. 


III. EXTENDED NIKIFOROV-UVAROV METHOD 


This method is improved for solving any second order differential equation which has at 
most four singular points. The improvement is achieved by changing the boundary conditions 
of NU method. An equation provided that the boundary conditions of extended NU method 
are satished, is reduced by an appropriate coordinate transformation into the following form; 


= 0 


CTAZ 


( 8 ) 


where Te{z), cre(z) and (Te(z) are polynomials of at most second, third and fourth-degrees 
respectively. By choosing a suitable function (j)e{z) and taking; 


tp(z) = 4>^(z)v{z) 


( 9 ) 


Eq.® is converted to the following form: 


y 


'U) + ( 


o 0e(^) ^ T-, 


(j)e{z) ae(z) 


y\z) + 


0e(^) I 0e(^) A(^) I ^ 

MZ) MZ) CTeiz) ^ ’ ■ 


( 10 ) 
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To obtain a simplified form of Eq. fllOp the coefficients of y{z) and y'{z) are arranged with 
newly defined polynomials. The coefficient of y'{z) is taken in the form Te{z)/ae{z)-, 


O 0e(^) I ^e{z) ^ Te[z) 

(peiz) (Je{z) ae{zy 


( 11 ) 


where Te{z) is a polynomial of at most second degree. The most regular form of this equation 


is obtained as follows; 


where 7ie{z), provided that 


0e(^) ^ T^ejz) 
y^z) ae(z) 

Te{z) = fe{z) + 2'Ke{z), 


( 12 ) 

(13) 


is a polynomial of at most second-degree. In addition for the coefficient of y{z), Eq. flTU]) is 
also rearranged as follows; 

</>e(^) I ^e{z) ^e(z) (fe(z) , . 

(j)e{z) (j)e{z)ae{z) a^(z) CT^(z)' 

where (fe(z) is a polynomial of at most fourth degree and it is defined as; 


o-yz) = o-e(z) + 7r^(z) + 7ryz)(Tyz) - o-'^z)) -f 7r'yz)ae(z). (15) 


After substituting the right-hand sides of Eq. llTTl) and Eq. (lT4)) into Eq. llTOl) the following 
equation can be obtained: 


/// ^ ^ , ^e[Z , , ^ 

If the polynomial (Je(z) is divisible by cre(z), i.e., 


(16) 


ae{z)lae{z) = h{z) (17) 

where h{z) is a polynomial of degree at most one, Eq. (lT^ is reduced to the following form; 

ae{z)y''{z) + Te{z)y'{z) + h{z)y{z) = 0. (18) 

If the definition ae{z)/ae{z) = h{z) is also used in Eq. flTSl) . the second order homogeneous 
equation is obtained as 

T^l{z) + Tle^Zyfe^z) - a'^z)) + d^e{z) - {h{z) - Tt'( 2))cTe(2) = 0, (19) 


where 

h{z)-'K'yz)= g{z). (20) 
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Using the roots of Eq. ffT^ . TVe^z) can be obtained as follows; 


vr. 


^ CTe(^) - Tejz) ^ 


,a'^{z) - Te{z) 


y -(Te{z) +g{z)ae{z). 


( 21 ) 


To determine all of possible solutions of the polynomial TTei^), ^^e polynomial g{z) under the 
square root sign must be known explicitly. Since 7Te{z) is a polynomial of at most second- 
degree, the expression under the square root sign must be square of a polynomial of at most 
second-degree. Therefore, the polynomials g{z) which satisfy this condition are chosen and 
after determination of g{z), the polynomial Tieiz) can be obtained from Eq. (l?Il) . Then Te{z) 
and h{z) can be found from Eq. flT^ and Eq. fl20|) respectively. 

In order to generalize the solutions of Eq. dTSj) . this equation is differentiated once and a 
compherensible form is obtained: 


ae{z)y^^\z) + {Te{z) + a'yz))y"{z) + {y{z) + h{z))y'{z) + h\z)y{z) = 0. (22) 


This equation is a third-order homogeneous differential equation with polynomial coefficients 
of degree not exceeding the corresponding order of differentiation. Since all derivatives of 
Eq. fl22l) have the same form, it can be differentiated n times by using the new representation 
|/W(^) = Vn{z)] 


((n -F l)h'{z) + 


ayz)vy\z) + {t^z) + {n + l)a'yz))v'^{z) + 

((n + l)r'(z) + ^^ <(^) + Kz))v'n{z) + 

n(n + l)_„^^^ , +1)(?^-1)_(3) 


-r'M + 


(Tl \z))Vn{z) = 0 


2 ^ 6 
When the coefficient of Vn{z) is equal to zero, polynomial hn{z) is dehned as 

, / N ’^ // \ ~ 1) /// \ ^ 

hn{z) = - ^ -CT.(z) + C„, 


(23) 


(24) 


where Cn is an integration constant. Therefore Eq. fl23ll is reduced into the following form: 

ff.(z)«®(z) + (t.(z) + (n+ 1 )(t;(z)^k;;( 2 ) + 

((^ + 1)<U) + d^<(z) + (z) = 0. (26) 

Eq. fl25ll has a particular solution of the form y{z) = yn{z) which is a polynomial of degree 
n. To obtain the eigenvalue solution, the relationship between h{z) and hn{z) must be set 
up by means of Eq. fl20|) and Eq. ( 
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IV. SOLUTION OF HEUN EQUATION BY EXTENDED NIKIFOROV-UVAROV 
METHOD 


Heun equation in which Gauss hypergeometric, confluent hypergeometric, Mathieu, Ince, 
Legendre, Laguerre, Bessel functions are involved, is a general second order linear differential 
equation: 

6 e dw a(5z-q _ 

^ + h + —+ ( 26 ) 

where 7 , 5, e, a and f3 are constant parameters which satisfy the Fuchsian condition; e = 
q; + /3 — 7 — (5 + 1. The complex number q is called accessory parameter. Eq. fl26p has four 
regular singular points at z = 0, z = 1, z = 00 and z = a js]. The most widely studied 
solutions of Heun’s equation literature are Heun polynomials which could be classified into 
eight classes. These solutions are written in the form; 


Hp{z) = z'^^i^z — iy^{z — aY^p{z) 


(27) 


where p{z) is a polynomial and (Ji, cr 2 , as are each one of the two possible exponents at the 
corresponding si^ularity. The classification of Heun polynomials is demonstrated in a table 
that is given as [9|; 


Class ai a 2 


0-3 


a 


/3 


I 

H 


0 0 

1-7 0 


0 

0 


—n 

7 — n — 1 


e + 7 + (5 + ? 7 , — 1 
n + 6 + e 


HI 

0 1 — (5 0 6 — n — 1 

e + 7 + n 

IV 

1 — 7 I — 5 0 (5 + 7 — n — 2 

e + n + 1 

V 

0 0 1 — e e — n — 1 

5 + 7 + n 

VI 

1—7 0 1 —e e+ 7 —n —2 

5 + n + 1 

VH 

0 1 — 51 — e S + € — n — 2 

7 + n + 1 

VHI 

1 — 7 I — 51 — ee + 7 + 5 — n — 3 

n + 2 


Table-I; Classification of Heun Polynomials 

Heun equation yields the boundary conditions of extended NU method. Comparing Eq. ( 
and Eq.dH]), the relationship between polynomials in the extended NU method and the 
parameters in Heun equation is set up as follows; 


Te{z) = ^{z — l){z — a) + dz(z — a) + ez(z — 1 ) 


(28) 
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ae(z) = z(z - l){z - a) 

cre(z) = {a[3z — q)z{z — l){z — a). 


(29) 

(30) 


Substituting these polynomials into Eq. fl^ 'n'eiz) polynomial can be obtained as; 


(1 - j){z - l){z - a) + (1 - 6)z{z - a) + (1 - e)z{z - 1) 

^e[z) = - - -± 

^{[(1 - 7){z - l){z - a) + (1 - 6)z{z - a) + (1 - e)z{z - l)f - 

4:{a(3z — q)z{z — l){z — a) + 4:g{z)z{z — 1)(2; — a)| 

Since 7r{z) is defined as a second order polynomial, the expression under the square root 
sign has to be the square of a second order polynomial. Hence all of the possible solutions 
of 71 (z) corresponding to certain values of polynomials g{z) can be listed as follows: 


TTeiiz) = (1 - j){z - l){z - a) + (1 - 6)z{z - a) + (1 - e)z{z - 1) 

T^e2{z) = 0 


for gi{z) = a(3z — q 


(31) 

(32) 


TTe^{z) = {1 -'y){z - l){z - a) (33) 

T^e^iz) = (1 — S)z{z — a) + (1 — e)z{z — 1) (34) 

for g 2 {z) = a^z - g - (1 - 7)[(1 - S){z - a) + (1 - e){z - 1)] 


Tiebiz) = (1 - 'y){z - l){z - a) + (1 - 6)z{z - a) (35) 

7re6(2;) = (1 - e)2;(2; - 1) (36) 

for g 3 (z) = a^z - q - (1 - e)[(l - 'y)(z - 1) + (1 - 6)z] 

7Te7{z) = (1 - 'y){z - l){z - o) + (1 - e)z{z - 1) (37) 

■Keniz) = (1 - 5)z{z - a) (38) 
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for gi{z) = afiz — g — (1 — 5)[(1 — ^)[z — a) + (1 — e)z\ 


Starting from vrei(z) in Eq. fl3ip the analytical solution can be written in details as follows: 
h{z) can be obtained from Eq. fl^UD : 

h{z) = a(3z - q + 2z[3 - (e + 7 + 5)] - (1 - 7 )(a + 1) - (1 - 6)a - (1 - e), (39) 

Then from Eq. (l24)) polynomial hn{z) is achieved as; 

hn{z) = -n[n + 5 - (e + 7 + 5)]z + 

Ti 2 

— [{a + 1)(2 — 7 ) + a(2 — 5) + (2 — e) + -{n — l)(a + 1)] + Cni, (40) 

By making use of the equality h{z) = hn{z) the eigenvalue equation and q accessory param¬ 
eter could be calculated and are given by; 


a(3 = {n + 2){e + ■y + 6 — n — 3) (41) 

Tt 2 

= 2 “ 7 ) + «(2 — <^) + (2 — e) -(n — l)(a 1)] -h 

(1 — 7 )(a -|- 1) -|- (1 — 6)a -|- (1 — e) -|- Cni- (42) 

The eigenfunction solution can be also acquired as follows; Firstly the function 0e(^) is 
calculated from Eq. flT^ : 

(t)e{z) = z^~'^{z-lf~\z-af~^. (43) 

Then a polynomial solution of degree n of Eq. fl2^ is obtained and demonstrated by Pi{z). 
Substituting these into Eq.(l9]), the eigenfunction can be achieved analytically; 

iJi{z) = (t)e{z)y{z) = z^-'^iz - 1)^“^(^ - aY~^pi{z). (44) 

This solution corresponds to the class VIII Heun polynomial j^. Since the above mentioned 
procedure is valid for all of the polynomials T^e{z) given in Eg. (132113^ in the eigenvalue and 
eigenfunction solutions, one can directly write down these solutions: 

For Eq. fl32D : 

a(3 = —n(e + 7 -|-( 5 -|-n — 1) 

Tl 2 

—9 — ■^[(^ + 1)7 + + € + —(n — l)(a + 1 )] + C,i2 
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( 46 ) 

( 46 ) 


- 1)°(2 - 


( 47 ) 


This solution corresponds to class I Heun polynomial [9|. 


For Eq. fl3^ : 


CK/S = (n + (5 + e)(7 — n — 1) (48) 

Tt 2 

—q = — [(a + 1)(2 — 7 ) + 0(5 + e + - (n — l)(a + 1)] + (1 — 7) ((3(5 + e) + (49) 


^{z - l)^{z - a)^p3{z) 

which corresponds to class II Heun polynomial j^. 

For Eq. flSTl) : 

CK/S = ((5 + e — n — 2)(7 + n + 1) 

Ti 2 

—q = —[{a + 1)7 + a (2 — (5) + (2 — e) + -{n — l)(a + 1)] + 

7 [a(l - (5) - (1 - e)] + 


(50) 


(51) 

(52) 


= z°{z - 1 )^"^(^ - ay-^p 4 iz) 
which is the corresponding class VII Heun polynomial jol- 


For Eq. (|35i) : 

a/5 = ((5 + 7 — n — 2)(e + n + 1) 

Tt 2 

-g — —[(a + 1)(2 — 7) + a(2 — 5) + e + -(n — l)(a + 1)] + 

(1 — 7)(a + e) + (1 — 5)a + (7„5 


(53) 


(54) 

(55) 


namely the corresponding class IV Heun polynomial 9|. 


(56) 
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For Eq. fl36D : 


a{3 = (5 + 7 + n)(e — n — 1) (57) 

Tl 2 

—q = ^[(o + 1)7 + (i5 + e — 2 + —(n — l)(a + 1)] — (2 — 7)(1 — e) + (58) 

= z^{z-l)^{z - ay~^Pfi{z) (59) 

Eq. fl5^ is the corresponding class V Heun polynomial {q]. 


For Eq. fl37|) : 


a(3 = (e + 7 — n — 2)((5 + n + 1 ) 

Tl 2 

—q = —[{a + 1)(2 — 7 ) + 0(5 + (2 — e) + -{n — l)(a + 1 )] + 

(1 — 7) (® “1“ (1 — "h Cn7 


(60) 

(61) 




This solntion is the corresponding class VI Henn polynomial 




For Eq. fl3^ : 


a(3 = (e + 7 + n)(5 — n — 1) 
71 2 

-9 = 21 (“ + I )'' + “(2 - f + 3 (" - 1 )(“ + 1)1 - 

(2 ~ 7)(1 ~ + C„s 


(62) 


(63) 

(64) 


tp 3 (z) = z"(z - 1 )‘ ‘(z - a)"ps(z) 


Eq. fl65l) gives the corresponding class III Heun polynomial 


a. 


(65) 


For the existence of polynomial solutions of Heun equation, one of the parameters a 
or f3 must be chosen conveniently. For instance, for the class VHI Heun polynomial these 
parameters are chosen tobea = e + 7 + (5 — n — 3 and /? = n + 2 j^. In the solution which 
is obtained by extended NU method, the parameters a and (3 can be determined exactly. 
Also the accessory parameter q is specihed exactly up to a constant of integration. 
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V. SOLUTION OF CONFLUENT HEUN EQUATION BY EXTENDED NIKIFOROV- 
UVAROV METHOD 


The ’confluent Heun equation’ is derived from the Heun equation by the coincidence of 


singu 

Ince 


arities at z = a and z = oo. The general confleunce process is discussed in details in 


in Ref. 


id • Also derivation of the confluent Heun equation (CHE) is given as an illustration 
d]. The most general form of CHE is expressed as; 




Tt ^ ^ - zsy 

" Bi A D 


2 r 


E 


E 

j= 

4 

E 


E. 


3j 


'3j 


y'{z) + 


'^{z) = 0 


( 66 ) 


where the location of singular points and coefficients are arbitrary. The simplest uniform 


shape of Eq. (l66D can be written as 


ill; 


y''{z) + (^a + 


y + 1 7 + l\ ^ ^ W ^ n 

T/u) + (t + 7-7)l/(U = 0- 


z - 1 


z z — 1 


(67) 


This equation has regular singularities at z = 0 and z = 1, and an irregular singularity at 


z = oo. The solution of Eq. 
by confluent Heun function 


ii in the vicinity of the regular singular point z = 0 is given 

3 


Hc{a, y, 7 , 6, r],z) = J2 7 , V, z)z'^, 


( 68 ) 


where 


n=0 


6 = fi + u-^{y+ 'J+ 2) 
r] = ^{y + 1) - /i - ^(/3 + 7 + /37) 


(69) 


and the coefficients Vn are given by three-term recurrence relation j^. This solution can be 
reduced to a confluent Heun polynomial of degree n provided that the following conditions 
are satisfied: 

fi + u = —na (70) 


together with the vanishing of the (n -|- l)x(n -|- 1) tridiagonal determinant given by 13| : 
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At -91 (1 + / 3 ) 


0 


0 


0 


0 


no fi — q 2 + a 2(2 + /3) . . . 0 0 

0 (n — l)o fj. — qs + 2a . . . 0 0 


0 

0 


= 0 . 


0 0 0 ... II — Qn-l + {n — 2}a (n — l)(n — 1 +/3) 0 

0 0 0 ... 2q 11 — Qn + {n — l)o n(n + /3) 

0 0 0 . . . 0 a II — qn+l + no 


The CHE which is given by Eq. flhTll . could analytically be solved by extended NU method 
as well. The solution process begins by the comparision of the boundary conditions of ex¬ 
tended NU method with the degrees of polynomial coefficients in CHE. After this compari¬ 
sion the polynomials in the extended NU method can be associated with the parameters in 
the CHE as follows: 


fe{z) = az{z - 1) + (/9 -F l){z - 1) -1- (7 + l)z (71) 

ae{z) = z(z - 1) (72) 

ae(z) = 1(/U + I 2 )z -/i]z(z - 1). (73) 

Subsituting these polynomials into Eq. (|^ . the polynomial 7re(z) can be obtained: 

-az(z - 1 ) - /3(z - 1 ) - 72 
= - 2 -^ 

- 1) - - 1) - izf - 4[(/i + u)z- g{z)\z{z - 1). (74) 


Since 7re{z) is a second degree polynomial, the expression under the square root must be a 
second degree polynomial. For this purpose polynomials g{z) are choosen as follows; 


vrei(^) 


for gi{z) = (/i -h i')z — fi 


—az{z — 1) — I3{z — 1) — 7 ^ 
T^e2{z) = 0 


(75) 

(76) 


= --iZ 

T^e^kA — —Q:^(^ “ 1) ~ ~ 1) 


13 


(77) 

(78) 






for g 2 iz) = {fi + u — ^a)z — /r — 7 /? 


ne 5 {z) =-/3{z - 1) --fz (79) 

T^esiz) = -az{z - 1) (80) 

for gsi^z) = {fj. + u — af3 — ^(y)z — /i — a(3 

Tferiz) = -/^(z - 1 ) 

TTesiz) = -az(z - 1) -'yz 

for g 4 (z) = (ju + ly — a/3)z — ju + a/3 — I3j 

For TTeiiz) in Eq. flTSl) . the analytical solntion is given as follows: polynomial h{z) can be 
obtained from Eq. fl20D : 

h{z) = + ly — 2a)z — /i + a — /3 — 7 . (83) 

polynomial hn(z) is achieved from Eq. fl2^ when the coefficient of Vn(z) is eqnal to zero as 
mentioned in Section.3. Bnt in the CHE problem a(z) is a second degree polynomial snch 
that a^^^(z) in the coefficient of Vn(z) is eqnal to zero. Then hn(z) is: 

Tl Ti 

hn{z) = ~y^T'{z) + Cn = — — (2az + 2 — q; + /3 + 7 ) + Cn- (84) 

By eqnating the right hand sides of Eq. flH^ and Eq. fIMll . eigenvalne eqnation and parameter 
/i (with a difference np to an integration constant) are achieved; 

fi + u — 2a = na (85) 

Tl 

—/i = + a + /3 — y') — a + /3 + 'y + Cni- (86) 

For the eigenfnnction solntion of CHE, (peiz) and the polynomial part of the solntion can 
be specihed from Eq. flT^ and Eq. fl2^ respectively. The fnnction (pe{z) is; 

= e-^^z-^{z-l)-\ (87) 


( 81 ) 

( 82 ) 
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If the polynomial solution is indicated by pi{z), by using Eq.([9]) the eigenfunction solution 
can be obtained completely; 


= (i)e{z)y{z) = e ^{z - 1 ) '^pi{z). ( 88 ) 


Since the solution process given for TTei{z) will be same for all of the polynomials 7ie{z) given 
in Eg. 07611821) ■ the results can be written down directly: 

For Eq. fl76D : 



/i + z/ = —na 

(89) 


Tl 

Ai=~(2 — Q: + /3 + 7)+ Cn2 

(90) 


i^2{z) =P2{z). 

(91) 

For Eq.fITTD: 


/i + n — ya = —na 

(92) 


71 

—p = — a + /I — 7 ) + 7(1 + /I) + Cn3 

(93) 


^z{z) = {z- 1 )"> 3 (^)- 

(94) 

For Eq.flTHl): 


p + 1 / — 'ya — 2a = na 

(95) 


Tl 

—/i = — "2 + 0 — P l) — Q' + /3 + 7/5 + C„4 

(96) 



(97) 

For Eq.(|7^: 


/i + z/ — 7 a — (3a = —na 

(98) 


Tl 

—p = + a — /5 + 7 ) — a + /5 + 7/5 + Cn5 

(99) 


= Z-^{z - l)"> 5 (z). 

(100) 

For Eq.fISOl): 


p + u — 7 a — /3a — 2a = na 

(101) 


Tl 

—p = — "2 + q: + /5 + 7) — Oi — al3 + CnQ 

(102) 


Mz) = e-^^pe{z). 

(103) 
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For Eq. flHTD : 


yU + z/ — jSa = na 

Ti 

—/i = + a — /3 + 7 ) — a(5 + / 3(7 + 1 ) + Cni 

^j(z) = z~^p7{z). 


For Eq. flH^ : 


jj, + u — f5a — 2a = na 

Ti 

—/i = -\-a + l3 — a — a(3 + 7(1 + /3) + C*,i 8 


(104) 

(105) 

(106) 


(107) 

(108) 
(109) 


VI. APPLICATIONS 

1. The Coulomb problem on a 3-sphere 

In Ref.[^, an analytical solution of the Schrodinger equation for the Coulomb problem on 
3-sphere in generalized parabolic coordinates is given. After using the method of seperation 
of variables for the wave function in the Schrodinger equation, they obtain the Heun equation 
in the following form; 


HI" + (1 + ^ + = 0 . 

z-l z+V z{z-l){z+l) 


( 110 ) 


where 


a = 1+ I m I 4 


F = 1 - \Jl + E + 

Afn dm I ^ I Al 

a/1 -1- E - Z7 - a/ 1 + 77 -f *7 


6 = 1-|- I m 


+ E - + E + 


( 111 ) 


In the above mentioned reference the eigenvalue solution of Eq. flllUp is obtained in the form; 


En = {n+ I m |)(n+ | m \ +2) — 


7 


4(n+ I m I +1)2 


( 112 ) 
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If the expressions given in Eg fillip are substitnted in Eq. fHsp and Eq. fHHP the same 
solution can be obtained by extended NU method; 

ab = —n{e + T + A + ?7, — 1) 

+ E - i'y + ^/l + E + i'y 


n , I I , , I 

(1+ I m I H-)(!+ I m 


-n{l — ^1 + E + *7 + 2 I m I +n + 2) 


En = {n+ I m |)(n+ | m \ +2) 


7^ 


4(n+ I m I +1)2 

(113) 


(1+ I m 


VI + -E - *7 - y'l + E + 


ab = (n + A + e)(r — n — 1) 

)(1+ I m I- - -) = 

(2 I m I +n + 2)(1 — + E + i'y — n — 1) 


En = {n+ I m |)(n+ | m \ +2) 


7 


4(n+ I m I +1)2 

(114) 


2. Two Coulombically repelling electrons on a sphere 


The problem of two electrons interacting via a Coulomb potential on a surface of a D- 
dimensional sphere of radius R is discussed in Ref. 1^ and the Heun equation is obtained 
for inter-electron wave function after introducing a dimensionless variable 2 ; = u/2R. ; 


(Pjj 

dz^ 


+ 


7 + 2 
iz 




z+1 


# -AR^Ez + 2R 
z — 1 \ dz z{z + l){z — 1) 


V = 0. 


It is also shown that this equation has polynomial solutions of degree n = 1, 2, 

n 


(115) 


(116) 


2 = 1 


stipulated that E and R take the values given by; 

n 


E = 


4R2 


(n + (5 - 1) 


fl=-5|2(«-l) + ^lE‘< 

^ i=l 


(117) 
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where ^ 1 , 2 : 2 , ■■■Zn can be determined from Bethe ansatz equations [iSj; 

r, ^ 1 


E' 


1(5-i) 1(5-i) 

_|_ Jl -j- ^ ^ 'y' 2 V ')> _ g 


Zi — Zj 


Zi 


Zi -h 1 


Zi - 1 


( 118 ) 


Furthermore they calculated the radius R and the energy E as R = E = 'y and 

R = |y^2(5 + 2) + E = for n = 1 and n = 2 states respectively. 

This problem could be analysed by an easier and faster approach when extended NU 
method is used. When the parameters 7 , 6, e, af3, q, a in Eq. fl2^ are replaced by I/ 7 ,1/2(5 — 
1 / 7 ), 1/2(5 — 1 / 7 ), —4:R?E, —2R, —1 respectively, one can obtain Eq. (11151) . After this match¬ 
ing, the energy E and the radius R (up to an integration constant) can be achieved by writing 
new parameters in Eq. flT5]) and Eq. fl46l) respectively: 


(119) 

( 120 ) 


- AR^E = -n(l /7 + 1/2(5 - I/ 7 ) + 1/2(5 - I/ 7 ) + n - 1) 

n 


E = 


4i?2 


(n -h 5 - 1) 

2R = a. 


For the state n = 1, the accessory parameter q = 2R must provide the determinant that 
admit polynomial solutions of Heun equation in Ref. jlT |: 
q -a'j 

= 0 . 

-ajS q + a{5 + 'y) + e + 'y 
For this problem this determinant has the following form: 

^ 0 . 

5 Cl 

From this determinant, the constant Ci can be obtained as follows; 


cl = Sh. 


( 121 ) 


Thus the radius R = 7 (by choosing the positive root since R is non-negative) and the 

energy E = 7 are obtained in agreement with those obtained in 15 1 . 

For n = 2 state analogously to state of n = 1, the accessory parameter q = 2R can be 

found from the solution of a determinant: 

q —07 0 

-a(3 q + a{6 +'y) + € +'y -2a{l + j) 

0 —(yjS — (7 T| ^ •if' 5) q 2(a -|- 1) -|- 2(a(5 + 7) T e + 7) 

which is given in Ref. ll| . For the new parameters, this determinant is transformed into the 

following form; 


=0 
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Cs 1/7 0 

2(1 + 5) C *2 2(1 + 7 ) “ 

0 2 + 5 C 2 

After solving for C 2 , we get; 


0 . 


^ 2 " = 2(5 + 2) + 


45 + 6 
7 


( 122 ) 


Therefore the radius R = + 2) + and the energy E = ^(-^(( 2 ')+ 2^+3 ^hich are again 


in full agreement with those obtained in 


ISj for the state n = 2. 


3. Hyperbolic double-well potential 

The solution of Schrodinger equation for a particle of mass m and energy E in & potential 
V {x) = is discussed in Ref.[^. One-dimensional Schrodinger equation for this 

particle is given by; 

+■)•(•'■ “■•'■Si)*""" 

where e = 2mE/h‘^, Uq = 2mVo/fi^ and s = x/d. In the aforesaid paper it is demonstrated 
that, this equation could be transformed into a confluent Heun equation; 

+ (« + 4^ + + (f + ^)i« = •> (124) 

where a = —d\/Uo, (3 = —idy/e, 7 = —1/2, /i = \{a{a + 2) + 2a(3 — — 1)), u = 

i(a+/3(/3+l)). To obtain the N-degree polynomial solutions of this confluent Heun equation, 
they have used the procedure which is reported in the hrst paragraph of Section 5 of this 
manuscript. After determining two termination conditions, the eigenvalue spectra for the 
symmetric (s) and antisymmetric (a) cases have been achieved as; 

£» =- 5 ^( 3 + 4Af-dvr;)' (126) 

e% = -^(5 + 4N-d^/uiy (126) 

This problem can be analysed by an easier and more efficient approach using the extended 
NU method. Since Eq. (ll24p has the same form with Eq. (1671) . the results which are obtained 
by extended NU method in Section 5, can be used directly for the eigenvalue spectra of this 
potential. Thus one can approach to the symmetric energy specta by using Eq. flH^ and 
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Eq. fll04p . When the parameters /i and v dehned for this problem are written in Eq. flH^ . 
one obtains; 

(X 

-(a + 3 + 2/3) =-na. (127) 

Subsituting a = —dy/Uo and (3 = —idy/e in Eq. fll27p the symmetric eigenvalue solution can 
be obtained exactly: 

1 / !—\2 

= + . (128) 

This solution is also achieved from Eq. fllOdp . The parameter /i which was obtained in Eq. fl^Up 
by extended NU method, can be modihed for this problem for the state iV = 1; 

/i = -(2 — CK + /3 + 7 ) + Ci 2 - (129) 

This result shows that the parameter /i can be determined up to an integration constant 
Ci 2 - Using the 2”'^ order determinant and the solution of related equation given in Section 
5, one could make a verihcation of this result. From this solution which is obtained via 
root-hnding methods, one could come to the conclusion: 

h = ^(2-a + /3 + 7)±^^(2-a + /3 + 7)2- 4a(l + (3). (130) 

When the integration constant C '12 is equal to {2 — a + (3 + 7 )^ — 4q:(1 + (3), our solu¬ 
tion will be the exact solution for parameter p. 

Proceeding as before the antisymmetric energy specta is achieved from Eq. fl92l) or Eq. fl98l) . 
Ater substituting parameters /i and v into Eq. fl9^ . it can be easily shown that 

OL 

— {a -|- 3 -|- 2(3) — 7 a = —na. (131) 

This allows us to hnd the following antisymmetric eigenvalue spectra: 

‘‘n ^. (132) 

The same result can be obtained upon solving Eq. (l98p as well. 


VII. CONCLUSION 

Nikiforov-Uvarov (NU) method which is developed for solving second order differential 
equations by special orthogonal functions, is an easy and elegant method. This method has 
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certain boundary comditions which are related to degrees of polynomial coefficients in the 
main equation of NU method. In turn, these boundary conditions bring some restrictions 
on the number of singular points of the relevant differential equations. Consequently, NU 
method will succeed only for equations which have at most three singular points. 

In this study, we propose that NU method is a method appicable for the equations which 
have at most four singularities. For this proposal the boundary conditions of NU method are 
changed and the extended NU method is obtained. Thereby, any second order differential 
equation which has at most four singularities, is reduced to the basic equation of extended 
NU method and analyzed systematically for eigenvalue and eigenfunction solutions. 

In this context, we have investigaed Heun and confluent Heun equations which yield the 
boundary conditions of the extended NU method and obtained exact eigenvalue solutions 
without using traditional power series methods. The eigenvalue solutions of three physical 
problems have been worked out in order to demonstrate the accuracy of the extended NU 
method and it has been observed that calculated energy eigenvalues are consistent with the 
former studies given in literature. 
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